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INTRODUCTION 

One  of  the  major  differences  between  steady-state  solution  techniques  for  the 
Navier-Stokes  equations  and  solution  procedures  for  either  inviscid  potential  or 
boundary  layer  problems  is  the  treatment  of  the  continuity  equation.  For  inviscid 
flow,  a  potential  function  is  determined  entirely  from  continuity.  The  pressure  is 
then  obtained  from  the  integrated  momentum  or  Bernoulli  equation.  For  boundary 
layers,  the  axial  momentum  and  continuity  equations  determine  the  velocities.  On  the 
other  hand,  typical  Navier-Stokes  solvers,  in  effect,  use  continuity  to  obtain  the 
density  (or  pressure),  and  the  velocities  result  solely  from  the  momentum  equations. 
For  large  Reynolds  number  steady  flows,  it  would  appear  that  such  a  procedure  is  in 
marked  conflict  with  both  the  asymptotic  inviscid  and  boundary  layer  theories. 

In  the  present  paper,  a  boundary  layer-relaxation  procedure  based  on  a  new 
composite-velocity  fohmilation  for  the  incompressible  Navier-Stokes  system  is  des¬ 
cribed.  The  equations  are  interpreted  and  numerically  approximated  to  reflect  the 
composite  nature  of  the  flow.  The  procedure  has  also  been  developed  independently 
for  subsonic  flow  [1] .  unlike  typical  Navier-Stokes  procedures  that  differ  signifi¬ 
cantly  from  their  incompressible  flow  counterparts,  the  present  developments  are 
essentially  identical  in  both  cases.  Moreover,  the  extension  to  transonic  flows 
should  be  direct. 

The  equations  are  written  in  a  body- fitted  orthogonal  coordinate  system  so  that 
arbitrary  geometries  can  be  treated.  Application  to  internal  and  external  flows  are 
discussed.  Specific  geometries  include  a  boat tail  simulator,  the  trailing  edge  of  a 
plate,  Joukowski  airfoils  and  a  curved  channel. 

In  its  final  form,  the  present  formulation  has  some  features  similar  to  the 
velocity-split  technique  due  to  Dodge ' (2] »  however,  this  resemblance  is  only  super¬ 
ficial.  In  the  present  analysis,  a  composite  representation  of  inviscid  and  viscous 
region  velocities  is  prescribed  in  the  spirit  of  matched  asymptotic  expansions.  The 
complete  Navier-Stokes  equations  are  solved.  No  simplifying  approximations  are 
required.  The  finite-difference  form  of  the  resulting  equations  are  solved  by  a 
coupled  strongly  implicit  procedure  described  previously  by  the  authors  13] . 

COMPOSITE  FORMULATION 

This  formulation  is  designed  for  the  calculation  of  large  Reynolds  number  pro¬ 
blems  with  a  dominant  flow  direction,  e.g.  the  (-direction.  The  gradients  are 


largest  in  surface  normal  (n)  direction.  The  flow  outside  the  thin  viscous  region  is 
essentially  inviscid  and  is  represented  by  a  potential  function  +;  therefore,  the 
following  composite  representation  of  the  velocity  field  to  reflect  the  matched 
asymptotic  boundary  layer* inviscid  behavior  is  prescribed. 

u  •  2-  (1  +  ♦.)  ■  u  0  ,  v  ■  ♦_  (1) 

n,  t  e  n,  % 


Substitution  of  these  expressions  into  the  ttavier-Stokas  equations  results,  after 
some  reorganization,  in  the  following  orthogonal  system  for  4,  a,  G: 
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For  internal  flows  G  (€)  is  determined  from  the  global  conservation  of  mass. 
In  these  equations 


u2*v2 


-  G(C> 


(5) 


and  is  similar  to  the  Bernoulli  or  total  pressure  for  the  inviscid  region.  However, 

G  is  not  assumed  constant,  but  is  determined  by  the  calculation  procedure.  In  an 
inviscid  irrotational  region,  0  ♦  1  and  continuity  (2)  reduces  to  the  well  known 
potential  flow  equation;  the  momentum  equations  (3,  4)  are  identically  satisfied  with 
G  *  constant  and  G  •  0.  For  internal  flows,  G  •  0  on  one  wall  and  G  is  determined 
by  mass  conservation.  In  the  viscous  region,  the  (  momentum  equation  determines  U, 
while  4  and,  therefore,  the  pressure  is  obtained  from  the  n  momentum  equation.  The 
two  basic  regions  pertinent  to  large  Reynolds  number  flows  are  appropriately  des¬ 
cribed  by  this  composite  set  of  equations.  This  method  of  defining  v  with  a 
'potential'  was  first  tested  for  the  flat  plate  boundary-layer  (U,4)  equations.  The 
solution  of  the  resulting  2-point  boundary  value  problem  reproduced  the  results  ob¬ 
tained  with  standard  methods  based  on  the  velocities  u  and  v.  It  should  be  pointed 
out  that  the  present  system  of  equations  can  be  used  for  the  solution  of  inviscid 
flown  provided  that  Be  ■  •  and  a  *  1  is  enforced  at  the  solid  surface.  G  •  0 
decouples  the  normal  moawntum  equation  from  the  axial  momentva  and  continuity  equa¬ 
tions  and  represents  the  interacting  boundary-layer  approximation  if  axial  diffusion 
is  neglected.  If,  in  addition,  ^  and  ^  in  equations  (2)  and  (3)  are  replaced 
with  their  respective  potential  flow  values,  the  usual  boundary-layer  approximation 

*V  dVOOTiTMe 


For  the  calculations  of  separated  flows  it  Is  essential  that  U  and  are 
coupled  in  order  to  eliminate  the  pressure  singularity.  G  which  represents  the 
Bernoulli  constant  in  the  in viscid  flow  can  be  treated  explicitly.  Finally,  the 
n -momentum  determines  G.  The  necessary  coupling  between  U  and  $  is  attained  with  a 
coupled  strongly  implicit  procedure  for  the  solution  of  the  algebraic  system. 

BOOHOARY  CONDITIONS 

In  the  present  formulation,  the  physical  boundary  conditions  are  represented 
directly  by  the  following  mathematical  boundary  conditions. 

At  a  solid  surfacei  o  ■  0>  <p  ■  0,  and  as  n  ♦  •,  $  **>  0,  u  1. 

n 

For  a  finite  body:  ♦  or  0  as  £  •*■  ±  •,  0(£,n)  «►  1  as  £  -•» 

and  0^(£,h)  *  0  as  £  ■*  •  . 

For  bodies  which  are  infinite  in  both  directions,  the  inflow  and  outflow  condi¬ 
tions  are  somewhat  ambiguous.  However,  these  conditions  must  be  appropriately  speci¬ 
fied  in  order  to  obtain  a  meaningful  solution.  For  the  boat tail  geometry  +  and  U 
are  prescribed  at  the  inflow,  while  0^  «►  0  and  ♦  0  have  been  applied  as  outflow 
conditions.  The  ♦  condition,  in  the  context  of  the  composite  formulation,  is  such 
as  to  eliminate  the  viscous -inviscid  interaction.  It  should  be  noted  that  these 
boundary  conditions  are  consistent  with  the  mathematical  character  of  the  equations 
governing  4  and  U.  For  example  no  slip  has  been  satisfied  through  the  boundary- 
layer  like  variable  0  and  not  through  4^. 

SOLOTIOH  PBOCBDORE 

The  governing  equations  have  been  discretized  using  second-order  accurate 
central-differencing  for  all  4  derivatives.  Central  or  boundary  layer-like  differ¬ 
encing  has  been  used  for  0  derivatives,  except  for  the  (hjU)^  term  in  the  continuity 
equation;  this  is  backward  differenced  throughout.  The  resulting  implicit  algebraic  | 
system  of  equations  has  been  solved  iteratively  using  a  coupled  strongly  implicit 
procedure  (CSX?) .  The  continuity  and  £ -momentum  equations  for  $  and  CJ  are  solved 
in  a  coupled  fashion,  while  the  n -momentum  equation  for  G  is  evaluated  iteratively. 

In  the  £  momentum  equation ,  G  is  treated  as  known  during  the  Iterations.  Although  G, j 
the  "inviscid”  total  pressure  is  evaluated  explicitly  from  the  n-momentum  equation, 

the  static  pressure  is  unknown  and  depends  upon  the  values  of  ♦_  and  *  .  Since  ♦ 

£  n 

Is  evaluated  implicitly  in  the  eocpled  algorithm,  this  implies  that  the  pressure  is 
also  treated  implicitly.  This  circumvents  the  separation  singularity.  Explicit 


artificial  viscosity  is  not  required  for  convergence  and  the  CSX?  allows  for  arbi¬ 
trarily  large  values  of  At  once  the  effects  of  the  inviscid  initial  conditions  have 
been  sufficiently  smoothed. 


COOP  USD  2x2  SOLUTION  ALGORITHM 

In  an  aarliar  paper,  rafaranca  [3],  tha  presant  authors  have  developed  a  coupled 
strongly  implicit  procedure  for  tha  stream  f unction-vortici ty  form  of  tha  Navier- 
Stokas  aquations.  This  algorithm  has  the  distinct  advantage  of  being  implicit  in 
both  tha  £  and  n  directions ,  as  wall  as  allowing  for  tha  coupling  of  all  tha  boundary 
conditions.  It  is  this  coupling,  that  eliminates  the  pressure  singularity  in  the 
(0,  ♦)  formulation.  Furthermore,  the  method  is  unconditionally  stable,  allows  for 
arbitrarily  large  At,  converges  faster  than  SOR,  LSOR,  ADI,  etc.,  and  is  relatively 
insensitive  to  grid  aspect  ratio.  The  discretized  version  of  the  equations  for 
(U,  ♦)  can  be  written  as: 

(A  +  P)  Vn+l  *  G  +  PV** ,  where  V  is  the  (0,  $)  solution  vector. 

P  is  chosen  such  that  (A+P)  can  be  decottposed  into  a  lower  and  upper  triangular  form 
having  a  sparsity  pattern  similar  to  the  original  matrix  A.  This  leads  to  a  solution 
algorithm  of  the  following  form: 


M-l 


i-l,j 


i-l,j 


where  n  is  the  iteration  index.  Although  the  coupling  accelerates  the  rate  of  con¬ 
vergence,  it  also  increases  the  storage  requirement  by  a  factor  of  two.  Considerable 
savings  in  storage  can  be  realized  by  re-evaluating  some  of  the  coefficients  (T.  .) 

j 

during  the  evaluation  of  U  and  4*  However,  this  is  achieved  at  moderate  additional 
computational  cost.  In  its  present  form  the  CSIP  is  slightly  different  from  the  one 
given  in  reference  [3] .  In  the  present  case  the  forward  and  backward  sweeps  have 
been  reversed  in  order  to  impart  a  certain  degree  of  marching  consistent  with  boundary 
layer  procedures.  As  detailed  in  reference  [3],  the  appropriate  recurrence  relation¬ 
ships  can  easily  be  obtained. 


RESULTS 

Laminar  flow  solutions  have  been  obtained  for  boattail  simulator,  Joukowski 
airfoil,  finite  plate  and  channel  geometries.  Reynolds  numbers  based  on  typical 
length  scales  ranged  from  103  to  105.  All  the  confutations  were  started  with  arbi¬ 
trary  initial  conditions.  For  the  first  40  iterations  At  •  1,  afterwhich  At  was 
increased  to  10*.  Explicit  artificial  viscosity  was  not  required  for  what  were 
effectively  steady-state  calculations. 

(a)  Boattail  Simulator:  A  typical  streamline  plot  for  an  axi symmetric  boattail 
geometry  is  shown  in  figure  1.  The  corner  angle  is  approximately  30*  and  the  body 


radius  varies  between  1  and  O.S.  With  1800  grid  points  (60  x  30,  12  on  the  boattail), 
and  Re-  7500,  convergence  to  10-3  is  achieved  in  about  125  iterations.  This  takes 
avout  10  minutes  on  the  Amdahl  470/V6  or  less  than  1  minute  on  the  Cray-1  computer . 
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(b)  Joukowskl  Airfoil:  X  variety  of  airfoil  thicknasaes  and  Reynolds  numbers 
have  been  considered.  For  low  Reynolds  numbers  and  small  thickness  ratio,  the  flow 
is  unsepar&tsd.  As  either  is  increased,  separation  regions  appear.  Typical  stream¬ 
line  pattern  and  separation  location  for  t/c  -  0.12,  0.1?  and  Re  ■  103  to  104  are  shown 
in  fig.'  2.  As  expected,  the  separation  point  moves  upstream  with  increasing  Reynolds 
number.  This  correlates  well  with  solutions  obtained  with  2nd  order  boundary-layer 
theory  [4] .  Small  changes  in  thickness  ratio  can  lead  to  large  variations  on  the 
recirculation  region.  The  solutions  are  oscillatory  for  large  t/c.  Re. 

(c)  finite  Flat  Plate:  The  flow  past  the  trailing  edge  of  a  finite  flat  plate 
has  been  extensively  investigated  by  triple-deck  interacting  boundary-layer  global 
relaxation  procedures  [5] .  In  figure  3  the  results  of  the  present  formulation  are 


compared  with  some  earlier  computations .  A  105  x  75  grid  has  been  adequate  to 
provide  reasonable  agreement. 

(d)  Internal  Flowt  Two  dimensional  straight  and  curved  channels  have  been  in¬ 
vestigated  by  the  present  technique.  The  curved  channel  was  generated  by  using  two 
streamlines  of  the  boattail  geometry.  A  variable  £  and  uni form- n  (60  x  30)  grid 
mbs  specified  for  this  calculation.  Uniform  inflow  conditions  are  prescribed.  The 
entrance  mass  flow  rate  was  0.4066  and  the  Reynolds  number  based  on  the  entrance 
channel  width  was  Re  *  2500.  Global  mass  conservation  was  insured  with  the  para¬ 
meter  G’  (£) .  Typical  wall  pressure  distributions  and  velocity  profiles  are  shown 
in  figure  4. 
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